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Braiding a flock: winding statistics of interacting flying spins 
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When animal groups move coherently in the form of a flock, their trajectories are not all parallel, 
the individuals exchange their position in the group. In this Letter we introduce a measure of 
this mixing dynamics, which we quantify as the winding of the braid formed from the particle 
trajectories. Building on a paradigmatic flocking model we numerically and theoretically explain 
the winding statistics, and show that it is predominantly set by the global twist of the trajectories 
as a consequence of a spontaneous symmetry breaking. 
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The collective behaviors observed in animal groups 
have attracted much attention in the biology, the math¬ 
ematics and the physics communities over the last 20 
years. Quantitative data analysis have established that 
the salient traits of collective motion are very well cap¬ 
tured by the dynamics of flying spins: persistent random 
walkers endowed with interactions akin to ferromagnetic 
couplings between their velocities [TH6]. This framework 
has been extensively exploited to rationalize structural 
and dynamical properties starting from the emergence 
of directed motion, to rapid (orientational) information 
transfer, see [B El El HE] and references therein. How¬ 
ever, beyond these spectacular results, the internal dy¬ 
namics of a flock, the relative motion of the individuals, 
remains scarcely investigated both experimentally and 
theoretically yet it is known to display non-trivial anoma¬ 
lous behaviors [HES]. 

In this Letter, we theoretically describe the mixing 
statistics of an archetypal flying-spin model. We first 
stress the intrinsic geometrical nature of the dynamics of 
particles in a flock and map this problem to the braiding 
statistics of their trajectories. We evidence the nontrivial 
statistics of the winding between pairs of motile-particle 
trajectories, which is a robust measure of their entangle¬ 
ment. This quantity displays spatial correlations at the 
population scale. We single out the reason for the non¬ 
trivial statistics and show that the spontaneous breaking 
of a rotational symmetry causes the global twist of the 
flock to chiefly rule the long-time winding fluctuations. 

Numerical flocking model. We build on a standard 
flocking model used to model compact groups akin to 
those observed in the wild m H H ED- N persistent 
random walkers, i = 1... TV, propel at a constant 

speed uo = 1. The dynamics of their orientation Pi(t) 
generically takes the form: 

^ = (I-P*Pi) • Fi({rj-,pj}j) (1) 

Rotational diffusion is accounted for by the uncorrelated 
Gaussian white noises of variance 2D. The interactions 
between self-propelled particles amounts to an effective 
torque which aligns the orientation p^ in the direction 


of Fi. The projection operator (I — PiPi) ensures that 
Pi lives on the unit sphere. We use a standard (metric) 
form of Fi to account for the flocking dynamics. Noting 
Tij = Ti — Tj it reads: 

= ZffA X P-J + ZffB X (2) 

* j&Ai » jgB, 

r is a relaxation time, which we henceforth set to r = 1. 
The first term in Eq. ([^ is a ferromagnetic term which 
promotes alignment with the mean direction of the 
neighbors lying in the sphere At of radius Ra = 1- The 
second term corresponds to attractive and repulsive in¬ 
teractions within the sphere Bi of radius Rb = 5, it is 
introduced to yield compact flocks m- Following mm, 
we assume that these interactions are attractive above 
a distance 2rc, and repulsive when r,. < Vij < 2rc. 
f{rij) = 1 - [rc/ir^j - Tc)]^ with Tc = 0.4. Eqs. Q 
and are solved numerically using an explicit Euler 
scheme. In all that follows, our only control parameter 
is the noise amplitude, D (see iiiiini for a comprehensive 
investigation of the phase behavior of this model). Below 
a critical amplitude ~ 0.3, the rotational symmetry 
of the particle orientation is spontaneously broken and 
collective motion emerges in the form of a flocking transi¬ 
tion (see the Supplementary Document [33]). A compact 
polar flock forms and moves along n(t) = {pi{f))i as 
exemplified in Fig. [^a) and in the supplementary video 
SI [34]. We restrict ourselves to this symmetry-broken 
phase and investigate the relative motion of the individ¬ 
uals within the flock. 

Quantifying the entanglement of the trajectories. 

We quantify the mixing dynamics of the flock by ex¬ 
ploiting a powerful toolbox that has been introduced in 
the context of Lagrangian mixing in fluids. The idea is to 
relate the mixing of an ensemble of moving particle to the 
entanglement of their wordlines [HHin]- A convenient 
measure of the entanglement of a bundle of trajectories is 
provided by a braid representation that we introduce be¬ 
low. Let us first define a convenient representation which 
disentangles the internal dynamics of the flock from the 
global turns of its mean direction of motion. We con- 
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FIG. 1: (a) Instantaneous positions and orientations of the particles in a compact polar flock (250 particles), and definition of 
the parallel-transported frame {fh, fjj, H). (b) World lines of 5 particles in the same polar flock. D — 2.6 x 10“^. (c) Definition 
of the crossing index, (d) Braid diagram associated to the world lines drawn in (b). 


sider an orthogonal frame {G{t), 772(1), n(i)) shown 

in Fig. Ba)- The origin is the center of mass G(t) of 
the flock. f)i{t = 0) is chosen arbitrarily; the f)^s are 
then parallel-transported along the trajectory of G in 
the course of the dynamics. If the flock were undergo¬ 
ing a rigid-body motion, the particle positions in this 
frame would be stationary. Conversely, mixing in the 
flock translates into the winding of the particle positions 
in {^ 1 ,^ 2 ) plane, see Fig. Bb) and supplementary movie 
S2 [34) . In order to quantify this winding dynamics, we 
observe the relative positions of the particles along a ref¬ 
erence axis, say rji. Particle winding is measured by 
observing the particle exchanges along f)i, and by as¬ 
signing an index e = ±1 to each crossing, as depicted in 
Figs. c) and supplementary movie S3 [34j. The 

sign of this index reflects the relative positions of the par¬ 
ticle in the orthogonal direction 7)2 as they cross along 7 )]^. 
More quantitatively, considering two particles i and j, we 
introduce their pair winding number, Wij{T), as the link¬ 
ing number between their world lines: Wij{T) = J 

where ej/ is the index of the crossing. The sum is per¬ 
formed over the crossings involving the particles i and j 
only, over a time interval T. This quantity has a clear 
meaning: it counts the number of turns of particle i 
around j (or, equivalently, of j around i) over a time 
T. The total winding in the flock between to and to + T 
is a measure of the entanglement of the world lines, and 
hence of the mixing: 


W{T) 


{id 


( 3 ) 


where N'p = N{N — l)/2 is the number of pairs. Im¬ 
portantly, W{T) does not depend on the distance be¬ 
tween the particles, only the index of the crossings and 
the times at which they occur matter. Therefore the 
world lines define a braid that can be drawn in the form 
of a normalized braid diagram. Fig. [^d) [H]. The to¬ 
tal winding corresponds to a topological invariant of the 
braid: W{T) = X)a where we now sum over 


all the crossings o this simplified representation. Prac¬ 
tically we choose the Artin representation of the braid 
word [B [B [2Tj . The braidlab library is used to 
compute both the pair and the total winding numbers. 
The topological nature of W (T) makes it a very robust 
measure of the flock mixing. We now carefully investigate 
its statistics. 

Winding statistics. 

The normalized distributions of W(T) are plotted in 
Fig-B a) for different values of the noise amplitude. For 
all the trajectory lengths, the total winding follows a 
Gaussian statistics with zero mean since the flock has no 
intrinsic chirality: clockwise and counter-clockwise wind¬ 
ings are equally probable. The winding distribution is 
fully characterized by its standard deviation (IF^(T))^, 
where the brackets denote the average over different ini¬ 
tial conditions. The winding fluctuations increase lin¬ 
early with the curvilinear length of the trajectories at 

short times: (IF^(T))^ oc T, and crosses over to a 

1 

diffusive regime where (IF^(T))^ oc Ts at long times. 
Fig. Bb). These first results would naively suggest a 
simple scenario. If the crossing events were uncorrelated 
the Gaussian statistics would readily stem from the cen¬ 
tral limit theorem as W(T) is the average of the crossing 
signs. In addition, the variance of W (T) would obviously 
grow in a diffusive manner. However, this appealing ex¬ 
planation is inconsistent with a deeper analysis of the 
data. Let us carefully study the correlations between the 
pair windings, which are quantified by: 


C^^(T) 
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{w,j{t)wki(T)) 

(wHT)) 


( 4 ) 

where {w‘^{t)) = A/p ^ is the variance of 

the pair windings. Given this definition, Cu,w(T) = 0 for 
uncorrelated WijS, and Gu,uj(F) = 1 when the pair wind¬ 
ings are fully correlated. Unexpectedly, we find that the 
correlation between the pair windings does not vanish at 
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FIG. 2: (a) Probability distribution of the total winding 

W(T), normalized by its standard deviation (T = 10^). The 
different noise amplitudes correspond to different polariza¬ 
tions of the flock 1^. Solid line; Gaussian distribution, 
(b) Standard deviation, (T))^^'^, as a function of time 
T, for different noise amplitudes, (c) Gorrelation function, 
defined by Eq. as a function of T. (d) Variance of the 
total winding normalized by the variance of the pair winding, 
plotted versus the particle number. Open circles: T = 10“'. 
Filled circles: T = 10®. Black squares: random walkers con¬ 
fined in a circular box of radius 20, diffusivity: 10. Time step: 
0.1. Dashed lines: slope —1. 


long times. Conversely, it increases and plateaus at a 
finite value, Fig. [^c), thereby ruling out the simple sce¬ 
nario sketched above. In order to check that this unex¬ 
pected behavior is not a finite-size artifact, we first note 
that Cu!w(T) ^ (IF^(T))/(rt;^(T)) in the large-A/’p limit 
and plot the ratio (W‘^{T))/{w‘^{T)) for flocks of differ¬ 
ent sizes iV, in Fig.j^d). Whereas short trajectories have 
indeed winding correlations C^w (T) that decay with the 
system size, as 1/N, the winding correlations of the long 
trajectories do not display any significant variations with 
N when increasing the particle number by a factor of ~ 8. 
In order to gain more insight into these tho opposite be¬ 
haviors, we computed the same quantity for the wordlines 
of independent 2D random walkers confined in a circular 
box. Cww follows the same 1/A^ nontrivial scaling ob¬ 
served for short flocking trajectories. We shall note that 
the finite time step of our numerical scheme regularizes 
the winding statistics of the random walkers and makes 
it possible to define its variance [23]. This second set 
of observations confirms that the saturation of C^w{T) 
at long time originates from extended correlations of the 
crossing events. 


A twist in the statistics. 

We now elucidate the physical mechanism responsible 
for the winding correlations. We first note that a global 
instantaneous rotation of the flock around 11 would re¬ 
sult in a fully correlated pair winding. We therefore sepa¬ 
rate the associated global trajectory twist from the total 
winding. Denoting by r]^ the position of particle i pro¬ 
jected in the observation plane the instanta¬ 

neous rotation rate of the flock is: 

(5) 

i 

Integrating over time, we define the global twist that is 
the number of turns of the flock around 11: 

^ rto+T 

Tw(T) = -y^ dtD(t). (6) 

We finally define the winding of the untwisted trajec¬ 
tories in the frame {G, r)i, 1 ) 2 ), obtained by rotating 
the parallel-transported frame (G, f);^, 7 ) 3 ) by an angle 
27r X Tw: 


W* = W- Tw. (7) 

Hence W* = where w*- = w^-Tw is the 

winding between particles i and j in this rotating frame. 
Within the braid picture, W* is found by factorizing out 
the global twist of the braid word [Ml [Mj . 

In Fig. I^a), we plot the standard deviation of W, Tw 
and W* versus T.The twist contribution dominates the 
total winding at long times and (W^) ~ (Tw^). This nu¬ 
merical fact is better understood by noting that the wind¬ 
ing in the rotating frame qualitatively follows the behav¬ 
ior displayed by confined random walkers. Above the re¬ 
laxation time r = 1 introduced in Eq. ([^, W*{T) follows 
the same diffusive evolution: {W*^{T)) oc T [Mj- More 
importantly, we also find the same 1 /N asymptotic scal¬ 
ing for Cw*w* ^ (T))/{w*'^(T)) showing that the 

spatial correlations of the displacements are short-ranged 
at all times, Fig. [^b). This result contrasts with the be¬ 
havior of the total winding in the parallel-transported 
frame where ^ (FF^(T))/(u>^(r)) hardly depends 
on N at long times. Fig. Hd)- We can therefore propose 
the following scenario: the long-range spatial correlations 
in the flock arise from the global rotation of the flock, 
hence from the global twist of the trajectories as clearly 
exemplified in the supplementary video S4 |34j . 

Hydrodynamic description of the trajectories ’ twist. 

The very origin of the global rotation of the flock roots 
from the spontaneous breaking of the rotational symme¬ 
try of particle velocities. This symmetry breaking gives 
rise to a soft orientational mode [27] which twists the tra¬ 
jectories in the parallel-transported frame at the entire- 
flock scale. We now lay out a more quantitative expla¬ 
nation by switching to an hydrodynamic description of 
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FIG. 3: (a) Standard deviations as a function of T: 

(open circles), (diamonds) and 

(Tw^(r))^^^ (filled circles), (VF^(T))^ ^ for confined random 
walkers (squares), same parameters as in Fig. 2. (b) Vari¬ 
ance of the total winding normalized by the variance of the 
pair winding computed in the twisting frame for different 
particle numbers. Open circles: T — 10“^. Filled circles: 
T — 10^. Black squares: confined random walkers. Dashed 
lines: slope —1. (c) Probability distribution of the twist nor¬ 
malized by its standard deviation {T — 10^). Solid line: 
Gaussian distribution, (d) The normalized twist diffusiv- 
ity: IlTw/(T^Tflock) decays quadratically with N. Solid line: 
slope —2. 


the flock viewed as a active-fluid drop. We use the con¬ 
ventional hydrodynamic framework first introduced phe¬ 
nomenologically by Toner and Tu m and later derived 
from microscopic theories HIMIIIS]. The fluid density 
and velocity fields are p{v,t) and v(r,t). We focus on 
strongly polarized flocks in which all particles follow the 
same average direction. The momentum equation lin¬ 
earized around the homogeneously polarized state takes 
the simple form [5J |3D] : 

ai(pv) -b A(fl • V)pv = -VP(p) -b rv2(pv) -b f, (8) 

where P is the local pressure, and f(r, t) is a Gaus¬ 
sian white noise with correlation (/^(r,t)/^(r', t')) = 
2DS{r — r')5{t — t')Sap- In this continuous limit, the 
global rotation rate, Eq. (§, is given by: 

n= ^ J d^r ^{r] X pv)-n. (9) 

Two comments are in order. Firstly, deep in the po¬ 
larized phase, the linearity of Eq. ([^ implies that the 


momentum fluctuations are Gaussian. After space and 
time integration, Eqs. © and ([^ imply that the twist 
also follows a normal distribution, in agreement with our 
numerical findings reported in Fig. [^c). Secondly, the 
damping of velocity fluctuations is set by the diffusive 
term rV^(pv), in Eq. Q (see pi 150] for more details). 
Hence the Fourier mode with wave-vector q decays in a 
time ~ (Tq^)”^. While small-wavelength perturbations 
are quickly damped, the large-scale fluctuations that oc¬ 
cur at the size of the flock, q ^ 1/Tflockj remain corre¬ 
lated over a time T ^ ^flock/^- This observation explains 
the time behavior of the trajectories’ twist fluctuations. 
At short times, T < T, the small-g fluctuations result 
in finite-time correlations in the rotation rate H. Gonse- 
quently, the twist fluctuations persist and undergo a “bal- 
istic” growth: (Tw^(r)) oc T^. At long times, T > T, 
all the Fourier modes have been relaxed, the correlations 
vanish, and one recovers the observed diffusive behavior 
for the trajectory twist: (Tw^(r)) ^ PtwE, Fig. (3)a. 
The scaling of the effective diffusivity with the size of 
the flock is computed in [3^: I?tw oc DT^LfiockN~‘^. 
This prediction agrees again with our numerical obser¬ 
vations, Fig. id). Now that we have elucidated the 
twist statistics, we finally explain why it dominates the 
winding in large flocks. Assuming that the flock density 
weakly depends on N, we find that (Tw^(T)) ~ 

As (FF*^)/(w*^) ^ N~^, Fig. ib), we conclude that the 
ratio between the local winding fluctuations and the twist 
vanishes as fV —>■ oo. 

In order to unambiguously prove that the soft rota¬ 
tional mode chiefly rules the winding statistics, we stress 
that all this phenomenology is lost in isotropic swarms. 
We show in the Supplementary Document |35j that when 
the noise amplitude is too strong to observe directed mo¬ 
tion, (VF^) deviates from (Tw^): increasing the noise re¬ 
sults in the decorrelation of the winding and the twist 
fluctuations [35| . Altogether, our numerical and analyti¬ 
cal results suggest a strong robustness of our main find¬ 
ings. The prominence of the twist fluctuations, leading to 
a coherent braiding of the flock trajectories, is expected 
to be qualitatively robust to the very details of the in¬ 
teractions. It solely relies on the spontaneous breaking 
of the rotational symmetry, and on its associated soft 
mode. It is therefore a direct consequence of the flocking 
transition and should be observed in all models yielding 
polarized flocks. 

We shall close with Letter from an experimental per¬ 
spective. In most of the situations, in the wild, external 
perturbations and fields explicitly break the rotational 
symmetry (e.g. gravity, predators, obstacles), or cause 
sudden collective turns as analyzed in [HEIHSa for star¬ 
ling flocks. The winding statistics should be an effective 
probe of the flock response to external bias without any 
a priori knowledge about the individual propulsion and 
interaction mechanisms. 
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Supplementary Information 


TRANSITION TO COLLECTIVE MOTION 


Upon decreasing the noise amplitude, the model defined by Eqs. (1) and (2) in the main text displays a transition 
to directed motion, upon decreasing the noise amplitude below a critical value ~ 0.2. The mean polarization. 
Ho = N~^ EiPil; increases from 0 in isotropic flocks to 1 in coherenty-moving groups, as shown in Fig. Both 
polar and isotropic flocks are compact: due to the attractive interactions, they do not span the entire simulation box 
but keep a finite size. A comprehensive characterization of the nature of the transition is provided in 
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FIG. 4: Mean polarization of the flock plotted versus the noise amplitude. 


LONG-TIME BEHAVIOR OF THE TWIST FLUCTUATIONS IN POLAR FLOCKS 


We detail the derivation of the scaling law for the twist diffusivity I?tw from the hydrodynamic description of the 
flock. For the sake of clarity, we introduce the momentum field V = pv. As we noted in the main text, its Fourier 
component Vq, with wave-vector q, is damped over a typical time ~ From Fq. (8), the time correlations 

therefore decay as: 

^ <5(q + q'), (10) 

where the indices a, P denote the spatial components of Vq. The smallest wave-vector being gmin 1/Aflock, the time 
correlations vanish is the limit where \t — t'\^ T AockA- regime, the momentum held is delta-correlated: 


OC “ t') ^a/3 (5(q + q')- (11) 

From this result, we deduce the autocorrelation of the global rotation rate. Fq. (9) provides the following expression: 

mW)) = ^jArJ (U^(t) Kit')), (12) 

where we recall that r/ = (I — IHI) • r and eapj is the fully antisymmetric Levi-Civita symbol. Moving to Fourier 
space and using Fq. 0, we have equivalently: 


D 


Ar / d^r' 




(U(t)fl(t')) 




6it-t'). 


(13) 







7 


readily provides the scaling 2 ?tw oc DT '^N for the effective diffusivity. Recalling that the persistence time of 

the twist is T we recover the scaling given in the main text: 

Ptw cx ^r'LflockAT-". (14) 


We therefore obtain (n(t)il(t')) ^ I^tw <5(t — t'). A simple dimensional estimate of the value of the integral in Eq. 13 


WINDING STATISTICS ACROSS THE TRANSITION TO COLLECTIVE MOTION 



D/D, 


D/D, 


FIG. 5: (a) Variance of the total winding (IF^(T) (blue open circles), and of the twist (Tw^(T) (red filled circles), plotted as 
a function of time T for a polarized flock {D = 0.02). (b) Same as (a) for an isotropic flock {D = 1.23). (c) Ratio between 
the variances of the twist and of the total winding at long times, plotted versus the noise amplitude, (d) Long-time correlation 
between the twist and the total winding, Cwtw = (lFTw)/[(lF^)(Tw^)]^'^^. 


In the main text, we restricted ourselves to polarized flocks. Here, we extend the analysis to higher noise amplitudes, 
above the transition to collective motion. As the flock has no well-defined direction of motion, we project the particle 
positions along a fixed direction (rather than in a parallel-transported frame). We compute the winding and the twist 
statistics as explained in the main text. We then compare, in Fig. the results obtained for polarized and isotropic 
flocks. We plot the variances (IV^(T)) and (Tw^(r)) in two extreme cases, see Figs, [^a) and (b). We also quantify 
the ratio (Tw^)/(1F^) at long times, as well as the correlations between the winding and the twist fluctuations, Cwtw, 
in Figs.j^c) and (d). Upon varying the noise across the transition to collective motion, the emergence of polar order 
is associated with a clear change in the winding statistics. 

• In polar flocks (at low noise amplitude, D < D, 0.2), we find that (Tw^) ^ (W^). The total winding and 
the twist are fully correlated, Cwt-w ~ 1- As they both follow a Gaussian distribution, these two quantities are 
statistically equivalent. 

• In isotropic flocks {D > D,), the total winding is well distinct from the twist as exemplified in Fig. [^b). In 
contrast with polar flocks, the ratio (Tw^)/(IU^) deviates from unity upon increasing the noise amplitude, see 













Fig. i c). This behavior is associated with a clear decay of the correlation between the winding and the twist 
in an isotropic flock, Fig. [^d). 

These results further confirm the scenario proposed in the main text to explain the winding statistics. In polar 
flocks, the total winding and the twist are both ruled by the same stochastic process: the spontaneous breaking of 
the rotational symmetry yields a spatially coherent rotation of the flock yet stochastic in time, i.e. a global braiding 
of the trajectories. This mechanism does not rely on the specific form of the interactions leading to polar order, it is 
therefore generic to all flocking models. By contrast, in isotropic flocks, this correlated rotational mode is suppressed as 
rotational symmetry is not broken. The residual winding and twist fluctuations solely arise from the weakly-correlated 
displacements of the individuals: in isotropic flocks, they correspond to different random processes. 
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